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Abstract 



We consider direct diphoton production in hadron collisions, and we compute the next- 
to-next-to-leading order (NNLO) QCD radiative corrections at the fully-differential 
level. Our calculation uses the qr subtraction formalism and it is implemented 
in a parton level Monte Carlo program. The program allows the user to apply 
arbitrary kinematical cuts on the final-state photons and the associated jet activity, 
and to compute the corresponding distributions in the form of bin histograms. We 
present selected numerical results related to Higgs boson searches at the LHC and 
corresponding results at the Tevatron. 
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Diphoton production is a relevant process in hadron collider physics. It is both a classical 
signal within the Standard Model (SM) and an important background for Higgs and new-physics 
searches. The origin of the Electroweak symmetry breaking is currently being investigated at the 
LHC by searching for the Higgs boson and eventually studying its properties. If the mass m# of 
the Higgs boson is low (m# < 140 GeV), the preferred search mode at the LHC involves Higgs 
boson production via gluon fusion followed by the rare decay into a pair of photons. Therefore, it 
is essential to count on an accurate theoretical description of the various kinematical distributions 
associated to the production of pairs of prompt photons with large invariant mass. Such task 
requires detailed computations of radiative corrections. 

In this paper we are interested in the process pp — > 77X (and the related process pp — > 77X), 
which, at the lowest order, occurs via the quark annihilation subprocess qq — > 77. The QCD 
corrections at the next-to-leading order (NLO) in the strong coupling «s involve the quark annihi- 
lation channel and a new partonic channel, via the subprocess qg — > 775. These corrections have 
been computed and implemented in the fully-differential Monte Carlo codes DIPHOX [lj, 2gammaMC 
[2] and MCFM [3j. A calculation that includes the effects of transverse- momentum resummation is 
implemented in RESBOS 0]. 

At the next-to- next-to-leading order (NNLO), the gg channel starts to contribute, and the 
large gluon-gluon luminosity makes this channel potentially sizeable. Part of the contribution 
from this channel, the so called box contribution, was computed long ago [5] and its size turns out 
to be comparable to the lowest-order result; for this reason, the box contribution is customarily 
included in all the NLO computations of diphoton production. The next-order gluonic corrections 
to the box contribution (which are part of the N 3 LO QCD corrections to diphoton production) 
were computed in Ref. [2] and found to have a moderate quantitative effect on the result of the 
'NLO+box' calculation. 

Besides their direct production from the hard subprocess, photons can also arise from frag- 
mentation subprocesses of QCD partons. The computation of fragmentation subprocesses requires 
(poorly known) non-perturbative information, in the form of parton fragmentation functions of 
the photon. The complete NLO single- and double-fragmentation contributions are implemented 
in DIPHOX pp. The effect of the fragmentation contributions is sizeably reduced by the photon 
isolation criteria that are necessarily applied in hadron collider experiments to suppress the very 
large irreducible background (e.g., photons that are faked by jets or produced by hadron decays). 
The standard cone isolation and the 'smooth' cone isolation proposed by Frixione |6J are two 
of these criteria. The standard cone isolation is easily implemented in experiments, but it only 
suppresses a fraction of the fragmentation contribution. The smooth cone isolation (formally) 
eliminates the entire fragmentation contribution, but its experimental implementation is still in 
progress. 

In this Letter we present the computation of the full NNLO QCD corrections to direct diphoton 
production in hadron collisions. We consider the inclusive hard-scattering reaction 

h + h 2 ^77 + X , (1) 

where the collision of the two hadrons, hi and I12, produces the diphoton system F = 77 with high 
invariant mass M 77 . The evaluation of the NNLO corrections to the process in Eq. ([1]) requires 
the knowledge of the corresponding partonic scattering amplitudes with X = 2 partons (at the 
tree level [7]), X — 1 parton (up to the one- loop level [8]) and no additional parton (up to the 
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two-loop level [9]) in the final state. The implementation of the separate scattering amplitudes in 
a complete NNLO (numerical) calculation is severely complicated by the presence of infrared (IR) 
divergences that occur at intermediate stages. The qx subtraction formalism [10J is a method that 
handles and cancels these unphysical IR divergences up to the NNLO. The formalism applies to 
generic hadron collision processes that involve hard-scattering production of a colourless high-mass 
system F. Within that framework [10J, the corresponding cross section is written as: 



d a (N)NLO 



^■(n)nlo ® da F Q + 



j -F+jets 7 CT 
aCT {N)LO ~ aa (N)LO 



(2) 



where daJ^^Q represents the cross section for the production of the system F plus jets at (N)LO 
accuracy!!], and da?^ LO is a (IR subtraction) counterterm whose explicit expression [12] is obtained 
from the resummation program of the logarithmically-enhanced contributions to qx distributions. 
The 'coefficient' H^ N ^ NLO) which also compensates for the subtraction of da^ LO , corresponds to 
the (N)NLO truncation of the pro cess- dependent perturbative function 
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The NLO calculation of da F requires the knowledge of 1-L F( ^\ and the NNLO calculation also 
requires 

The general structure of 'H F ^ is explicitly known [T5] : 1-L F ^ is directly obtained from the 
process-dependent scattering amplitudes by using a process-independent relation. Exploiting the 
explicit results of T-i F ^ for Higgs [TQl [14] and vector boson [15] production, we have generalized 
the process-independent relation of Ref. p3] to the calculation of the NNLO coefficient T-i F< ^ (this 
general result is presented in a forthcoming paper). Using this relation and the relevant scattering 
amplitudes [Hill [9], we have explicitly determined H F ^ for diphoton production. 

We have performed our fully-differential NNLO calculation of diphoton production according 
to Eq. (121). The NNLO computation is encoded in a parton level Monte Carlo program, in which 
we can implement arbitrary IR safe cuts on the final-state photons and the associated jet activity. 
The present formulation of the q^ subtraction formalism [TU] is restricted to the production of 
colourless systems F and, hence, it does not treat parton fragmentation subprocesses (here F 
includes one or two coloured partons that fragment). Therefore, we concentrate on the direct 
production of diphotons, and we rely on the smooth cone isolation criterion [6J. Considering a 
cone of radius r = \J (A77) 2 + (A0) 2 around each photon, we require that the total amount of 
hadronic (partonic) transverse energy Et inside the cone is smaller than ETmaxir), 



& 



Tmax{r) = e-fPj, 



1 



cosr 



1 — cos R 



(4) 



where p T is the photon transverse momentum; the isolation criterion Et < ETmax{ r ) has to be 
fulfilled for all cones with r < R. The isolation parameters are set to the values e 7 = 0.5, n = 1 
and R = 0.4 in all the numerical results presented in this Letter. We use the Martin-Stir ling- 
Thorne-Watt (MSTW) 2008 [16] sets of parton distributions, with densities and as evaluated at 
each corresponding order (i.e., we use (n+ l)-loop as at N n LO, with n = 0,1, 2), and we consider 
Nf = 5 massless quarks/antiquarks and gluons in the initial state. The default renormalization 



tin the case of diphoton production, the NLO calculation of dajj^o was performed in Ref. [IT] 
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(fin) and factorization (/ip) scales are set to the value of the invariant mass of the diphoton system, 
fJ'R — fJ'F — M 77 . The QED coupling constant a is fixed to a = 1/137. 

We apply typical kinematical cuts [17] that are used by the ATLAS and CMS Collaborations 
in their Higgs search studies. We require the harder photon to have a transverse momentum 
^harder > Q e y ? ^jjg f or the softer photon we demand p|? fter > 25 GeV. The rapidity of both 
photons is restricted to \y 7 \ < 2.5, and the invariant mass of the diphoton system is constrained 
to lie in the range 20 GeV < M 77 < 250 GeV. 



a (fb) 


LO 


NLO 


NNLO 


f>F = Hr = M 77 /2 


5045 ± 1 


26581 ± 23 


45588 ± 97 




5712 ±2 


26402 ± 25 


43315 ± 54 


y>F = Hr = 2M 77 


6319 ±2 


26045 ± 24 


41794 ± 77 



Table 1: Cross sections for pp — > 77 + X at the LHC (yfs = 14 TeV). The applied cuts are 
described in the text. 
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Figure 1: Invariant mass distribution of the photon pair at the LHC (\fs = 14 TeV): LO (dots), 
NLO (dashes) and NNLO (solid) results. We also present the results of the box and NLO+box 
contributions. The inset plot shows the corresponding K-factors. 

We start the presentation of our results by considering diphoton production at the LHC (y/s = 
14 TeV). In Table HJ we report the results of the accepted cross section at LO, NLO and NNLO. 
We have fixed fip = hr = jj, and we have considered three values of ///M 77 (/i/M 77 = 1/2, 1,2). 
The numerical errors estimate the statistical uncertainty of the Monte Carlo integration. 

We note that the value of the cross section remarkably increases with the perturbative order 
of the calculation. This increase is mostly due to the use of very asymmetric (unbalanced) cuts 
on the photon transverse momenta. At the LO, kinematics implies that the two photons are 
produced with equal transverse momentum and, thus, both photons should have Pj, > 40 GeV. 



3 



At higher orders, the final-state radiation of additional partons opens a new region of the phase 
space, where 40 GeV > p^? fter > 25 GeV. Since photons can copiously be produced with small 
transverse momentum (see also Fig. [2] and the related discussion), the cross section receives a 
sizeable contribution from the enlarged phase space region. This effect is further enhanced by 
the opening of a new large-luminosity partonic channel at each subsequent perturbative order. 
For example, at NLO the qg channel accounts for about 80% of the increase of the cross section. 
Therefore, it is not unexpected that a naive analysis of scale dependence (as presented in Tabled]) 
underestimates the size of the higher-order corrections. 
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Figure 2: Transverse-momentum distribution of the harder (left) and softer (right) photon at the 
LHC (yfi= 14 TeVj. 

In Figure [U we compare the LO, NLO and NNLO invariant mass distributions at the default 
scales. We also plot the gluonic box contribution (computed with NNLO parton distributions) and 
its sum with the full NLO result. The inset plot shows the K-factors defined by the ratio of the 
cross sections at two subsequent perturbative orders. We note that k nnl °/ nl ° [ s sensibly smaller 
than K iVLO//LO , and this fact indicates an improvement in the convergence of the perturbative 
expansion. In particular, the impact of the full NNLO corrections turns out to be reasonably 
moderate, with a K-factor, defined as the ratio between the NNLO and NLO+box distributions, 
of about K ~ 1.35. We find that about 30% of the NNLO corrections is due to the gg channel 
(the box contribution is responsible for more than half of it), while almost 60% still arises from 
the next-order corrections to the qg channel. The NNLO calculation includes the perturbative 
corrections from the entire phase space region (in particular, the next-order correction to the 
dominant qg channel) and the contributions from all possible partonic channels (in particular, a 
fully-consistent treatment of the box contribution to the gg channel . Owing to these reasons, the 
NNLO result can be considered a reliable estimate of direct diphoton production, although further 
studies (including independent variations of [Mr and /xp, and analyses of kinematical distributions) 
are necessary to quantify the NNLO theoretical uncertainty. 

In Fig. H] we show results on more exclusive observables: the pt distributions of the harder 
(left-hand plot) and softer (right-hand plot) photon. The statistical errors of the Monte Carlo 

■'"The calculation [5] of the next-order gluonic corrections to the box contribution indicates an increase of the 
NNLO result by less than 10% if M 77 > 100 GeV. 
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integration are at the per-cent level and hardly visible in Fig. [2J As previously anticipated, in 
the right-hand plot we observe that significant NLO and NNLO contributions to the cross section 
originate from the phase space region (25 GeV < p^? fter < 40 GeV) that, at LO, is kinematically 
forbidden by the asymmetric transverse-momentum cuts. In this low-pr region, the production 
mechanism of the softer photon is dynamically enhanced (the production probability is roughly 
proportional to a s x (p£r dor /p*? fter ) ln(^ arde 7p^ ofter ), if ^ ofter /^ arder ~ (^ ofter /M 77 ) 2 < 1), and 
this is responsible for a substantial part of the large higher-order corrections observed in Table HJ 
Fig. [1] and also in the px distribution of the harder photon. Comparing the px distributions of the 
harder and softer photon in the high-pr region (px ^ 50 GeV), we also observe that the distribution 
of the softer photon receives higher-order corrections that are sensibly smaller. This decrease is 
expected: if both photons have high px, the effect of the very asymmetric transverse-momentum 
cuts is reduced. 

We also comment on the Pt distribution of the softer photon in the region around the LO 
threshold (p|? ~ 40 GeV). Here the LO result has a step-like behaviour, and this necessarily 
produces [18] integrable logarithmic singularities at each subsequent perturbative order. The 
peak of the NLO and NNLO distributions at p^? fter ~ 40 GeV is an artifact of these perturbative 
instabilities. The instability is cured by all-order perturbative resummation, which eventually 
leads to a smooth pt distribution with a shoulder-like behaviour [18] in the vicinity of the LO 
threshold. This physical behaviour can be approximated (mimicked) in the NLO and NNLO 
calculations by smearing the distribution over a bin (with a sufficiently-large size) centered around 

^softer _ 4 Q Ge y 




Figure 3: Same as Fig. Ui but for Tevatron kinematics (as described in the text). 

In Figure El we present the invariant mass distribution for diphoton production at the Tevatron 
(a/s = 1.96 TeV). We require the harder and softer photon to have a transverse momentum larger 
than 17 GeV and 15 GeV, respectively. The rapidity of both photons is restricted to \y 7 \ < 1. We 
note that the increase from the LO to the NLO result is considerably smaller than in Fig. [TJ this 
is mostly due to the use of photon transverse-momentum cuts that are only slightly asymmetric. 
In the region where M 77 > 80 GeV, the relative impact of the box contribution is smaller than at 
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the LHC: this is a consequence of the higher values of parton momentum fractions, x, that are 
probed by Tevatron kinematics. Nevertheless, the NNLO corrections (which are dominated by the 
next-order correction to the qg channel) still increase the result at the previous order by roughly 
30%. 

We have presented the calculation of the cross section for diphoton production up to the 
complete NNLO in QCD perturbation theory. At the NNLO, all the contributions from the gg 
channel are included in a fully-consistent (and unambiguous) manner. Considering the illustrative 
isolation and kinematical cuts implemented in this paper, we find increasing effects of about 
30%— 40% with respect to computations at the previous perturbative order. Our calculation is 
directly implemented in a parton level Monte Carlo program. This feature makes it particularly 
suitable for practical applications to the computation of distributions in the form of bin histograms. 
A public version of our program will become available at a later stage. 
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